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Abstract 

We present analytic results for first order radiative QCD corrections to e + e~ 
annihilation processes into quarks and gluons under the special aspect of polarized 
final state particles. This aspect involves single-spin polarization, gluon polarization, 
and the correlation between quark spins for massive final state quarks. 



Contents 



1 Introduction 



2 The joint quark-antiquark density matrix 




2.1 The differential cross section 





2.2 The hadron tensoi 




2.3 The electro- weak form factors 


2.4 The polar angle dependence 





2 

3 
4 
5 
6 

7 

7 

8 
9 

3.4 Finite gluon energy cut| 12 



Unpolarized and polarized structure functions 

371 Born term contributions 



3.2 First order loop contribution? 



ITS" First order tree graph contributions 



Longitudinal spin correlation! 12 

171 Result for the correlation! 13 



4.2 Anomalous spin- flip contributions 13 



5 Gluon polarization 



14 

15 

5.2 Results for the gluon polarization| 15 



O The differential cross section 



6 Conclusion 



16 



Invited talk given at the Fiiiisika Instituut, 
Tartu Ulikool, Estonia, May 14th, 2002 



1 



1 Introduction 



Quark mass effects pfay an important roie in the production of quarks and gluons in e + e~ 
annihilations. Jet definition schemes, event shape variables, heavy flavour momentum 
correlations and the polarization of quarks and gluons are affected by the presence of 
quark masses for charm and bottom quarks even when they are produced at the scale of 
the Z° mass E3, BJ. A careful investigation of quark mass effects in e + e~ annihilations 
may even lead to an alternative determination of the quark mass values [|l], 0, [3], Q . 

There is an obvious interest in quark mass effects for ti production where quark mass 
effects cannot be neglected in the envisaged range of energies to be covered by the Next 
Linear Collider (NLC). The longitudinal polarization of massive quarks affects the shape of 
the energy spectrum of their secondary decay leptons, the longitudinal spin-spin correlation 
effects in pair produced quarks and antiquarks will lead to correlation effects of the energy 
spectra of their secondary decay leptons and antileptons. Quark mass effects are important 
in the m — > calculation of radiative corrections to quark polarization variables because 
residual mass effects change the naive no-flip pattern of the m = polarization results ||, [| . 

In this talk I present 0(a s ) radiative corrections to different polarization observables, 
namely the longitudinal and transverse single-spin polarization of the quark |7], [8], |||, 10 



the polarization of the gluon ITT], [TI|], and the longitudinal spin-spin correlation |13| , |I4 | 



for massive quark pairs produced in e + e~ annihilations, always including the polar angle 
dependence. As a byproduct of our calculations I discuss the m — ► limit and the role 
of the 0(a s ) residual mass effects. I explain how residual mass effects contribute to the 
various spin-flip and no-flip terms in the m — > limit for each of the three structure 
functions that describe the polar angle dependence. 



2 The joint quark-antiquark density matrix 

Let us start with the joint quark-antiquark density matrix d<sr a = da" x x , x , where Ai 
and A2 denote the helicities of the quark and antiquark, respectively. In the applications 
covered by this talk, only the diagonal case Ai = X[ and A 2 = X' 2 is of importance. The label 
a specifies the polarization of the initial 7*, Z or interference contributions thereof. This 
index will be specified later on by the labels U (unpolarized transverse), L (longitudinal), 
F (forward-backward), and the two labels A and / for the case of transverse polarization 
of the final state quarks. 

The diagonal part of the differential joint density matrix can be represented in terms of 
its components along the products of the unit matrix and the z-components of the Pauli 
matrix <tr 3 (tr 3 = p 1 <sr for the quark and <tr 3 = p 2 <ST for the antiquark, pi = Pi/\f>i\). One has 

da a = - (da a l ® 1 + da^a 3 ® 1 + da^H ® «r 3 + da^a 3 ® <r 3 ) (1) 
as formulated for longitudinal polarization, the first and the second Pauli matrices stand 
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for the quark and the antiquark, respectively. An alternative but equivalent representation 
of the longitudinal spin contributions can be written down in terms of the longitudinal 
spin components s{ = 2Ai and s l 2 = 2X 2 with s{,s 2 = ±1 (or s\,s 2 £ {Til})- The spin 
dependent parts of the differential cross section are given by 



da a (s{, s £ 2 ) 

Eq. (0) is easily inverted, 

da a 

da^ 
da^ 



- (da a + da^s[ + da^si + da^s[s e 2 ) . 



d<r a (TT) + da a (U) + <MiT) + da a (U), 
^q(TT) + da a (U) ~ rfo"a(IT) - da a (U), 
d<7a(tt) - da a (U) + da a (^) - dcr a (ll), 
da a (1^) - da a (U) - da a ([\) + da a {[[). 



(2) 



(3) 



The first line is the unpolarized contribution, the second line the contribution due to a 
single polarized quark, while the last line is the correlation between quark and antiquark. 
While all these contributions will be considered in the following, the polarisation of a 
single antiquark (third line) is not considered because it is equivalent to the case of a single 
polarized quark. 



2.1 The differential cross section 

According to Fermi's "golden rule" the differential cross section is given by 

da = 2n\T fi \ 2 dPS. 
Using conventional Feynman rules, especially 

-ieQfYn and - ieiyfip - a/T/ZTs) 
for the vertex of the fermion line with a photon or a Z boson, respectively, and 



-w 



-w 



and xz{q)- 

for the boson propagators with momentum q, one can write 



T fi = v{p 2 ) - ieQ f Y u{p 1 )—v(p + ) - ieQ^ «(p_ 



+v(p 2 ) [ - ie(v f Y - a f ^^)\u{p 1 )xz{q)-^v{p + ) 



u(p. 



(4) 



(5) 



(6) 



(7) 



where p_ and p_ are the moments of the electron and positron, pi and p 2 are the moments 
of the quark and antiquark, and q = P- +p+ = Pi +P2 is the moment of the boson. Possible 
polarisation degrees of freedom enter the squared matrix element by replacing for instance 
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u(pi)u(pi) by (j^Pi + m)(l + •y 5 'j u Si)/2 instead of summing over the polarisations of the 
final state, m is the mass of the quark (and antiquark). The squared matrix element reads 

\T fl \ 2 = g^H^ (8) 

where L l contains all elements due to electron and positron (lepton tensor), H^ v contains 
all elements due to the quark and antiquark (hadron tensor), and gij contains the elements 
due to the intemediate bosons (photon and Z boson) which mix up. These parts will be 
considered in turn. 



2.2 The hadron tensor 

The hadron tensor is determined by the hadron dynamics, i.e. by the current-induced pro- 
duction of a quark-antiquark pair which, in the 0(a s ) case, is followed by gluon emission. 
In the 0(a s ) case one also has to add the one-loop contribution. The index i = 1,2,3,4 
in Eq. (|S]) specifies the current composition in terms of the parity-even (for i = 1,2) and 
parity-odd {i = 3, 4) products of the vector and the axial vector currents according to 

rrl __ _(rrVV , ttAA\ tt2 _ \(ttVV _ ttAA\ 

In order to be specific, the unpolarized hadron tensor is e.g. given by the matrix 

H™=(qq\jV\0)(0\(JW\qq). (10) 
In the two-body case e + e~ — > qq the components of the cross section are given by 

„2„ 



r ] = ^ H « {P) (with v = y/l - 4mVg 2 ) (11) 



where (P) stands for unpolarized, single-spin or spin-spin correlation contributions. In the 
three-body case e + e~ — > qqg the hadron tensor components H£ p '(y,z) are related to the 
components of the differential cross section by 

S=iHi^^-4 (12) 

The two energy-type variables y = 1 — 2piq/q 2 and z = 1 — 2p2q/q 2 are used as kinematic 
variables. Note that the three-body helicity structure functions H^ p \y, z) have a different 
dimension than their two-body counterparts in Eq. (|T2"| ) which we indicate by explicitly 
referring to the (y, ^-dependence of the three-body structure functions. An example for 
the phase space in terms of y and z is shown in Fig. f|. 
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2.3 The electro- weak form factors 



The second building block = 1, 2, 3, 4) specifies the electro-weak model dependence 

of the e + e~ cross section. For the present discussion we need the components gn, g\2, #43 
and (744. They are given by 



9u 


= Ql - 2QfV e v f Rexz + {vl + a 2 e )(v 2 f + al] 


\Xz\ 2 


9 12 


= Q 2 f -2Q f v e VfRexz + {v 2 + a 2 e )(v 2 f - a 2 f 


)\Xz\ 2 


913 


= -2Q f v e a f lmxz, 




9i4 


= 2Q f v e a f Rexz-2(v 2 + a 2 e )v f af\xz\ 2 , 




921 


= q 2 f - 2Qtv e VfRexz + (v 2 - a 2 e )(v 2 f + a 2 f ) 


Xz\ 2 , 


922 


= q) -2Q f v e v f Rexz + {v 2 e - a 2 e )(v 2 f - a 2 f ) 


Xz\ 2 , 


923 


= —2QfV e af Im xz, 




924 


= 2QfV e afRexz-2(v 2 -a 2 e )v f af\xz\ 2 , 




931 


= -2Q f a e v f lmxz, 




932 


= -2Q f a e v f lmxz, 




933 


= 2Q f a e a f Rexz, 




934 


= 2Q f a e a f lmxz, 




941 


= 2Q f a e v f Rexz -2v e a e (v 2 + a 2 )\xz\ 2 , 




942 


= 2Q f a e v f Rexz - 2v e a e (v 2 f - a 2 f )\xz\ 2 , 




943 


= 2Q f a e a f lmxz, 




944 


= -2Q f a e a f Rexz + 4v e a e v f a f \xz\ 2 





(13) 



(14) 



(15) 



(16) 

where Xz(q 2 ) = gM%q 2 / \q 2 - M% + iM z T z ) , with M z and Y z the mass and width of the Z° 
(M z = 91.1887 GeV and V z = 2.487 GeV PJ) and g = G F {%^2ixa)- 1 « 4.49x lO^GeV -2 . 



Qj are the charges of the final state quarks to which the electro- weak currents directly 
couple; v e and a e , Vf and aj are the electro-weak vector and axial vector coupling constants. 
For example, in the Weinberg-Salam model, one has v e = —1 + 4sm 2 9w, a e = —1 for 
leptons, Vf — 1 — I sin 2 9w, a/ = 1 for up-type quarks (Qf = |), and Vf = — 1 + | sin 2 9w, 
cif = —1 for down- type quarks (Q/ = —3). The left- and right-handed coupling constants 
are then given by gi = v + a and gn = v — a, respectively. In the purely electromagnetic 
case one has gn = gy± = ^21 = #22 = Q 2 and all other g r i r = 0. The terms linear in 
Rexz and Imxz come from 7 — Z° interference, whereas the terms proportional to \xz\ 2 
originate from Z-exchange. 
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The generalization of the case where one starts with longitudinally polarized beams is 
straightforward and amounts to the replacement 



g u -> (1 - h h + )g u + (h - h + )g 4i g 4i -> (1 - h h + )g 4i + (h -h + )g u (17) 

where h~ and h + (—1 < < +1) denote the longitudinal polarization of the electron and 
the positron beam, respectively. Clearly there is no interaction between the beams when 
h+ = h- = ±l. 

2.4 The polar angle dependence 

As mentioned in the introduction, the polar angle dependence will be determined. The 
rest frame of the boson is a natural choice for the considerations. However, two coordinate 
systems are in use here, the event frame for the outgoing particles and the beam frame 
for the incoming particles. Transforming from one to the other system brings in a relative 
polar angle 9 (and in case of transverse spin vectors also a azimuthal angle (p which is 
averaged out normally). The polar angle 9 is given by the angle between the colliding 
electron beam direction and the outgoing quark direction. Transformed to the event frame, 
the lepton tensor L l leads to the different angular distributions. While for massless 
leptons the components L 2 and I? vanish, one can decompose the remaining lepton 
tensor components according to 

L 1 = *j + cos 2 9)n v + sin 2 9 II L +sin0cos0i7/j , 
a 2 

L 4 = ^-{cos9n F + sm9II A }. (18) 

Remark that 77/ and H \ depend linearly on siny? and cosy?. The matrices ilj/, iJ^, 
III , n F and n a are called projectors because in contracting the lepton tensor with the 
hadron tensor they project out the contribution of the hadron tensor to the different angle 
dependences. Explicit projectors are 

nl iv _ ? au VlVl jriiv _ VlVl jjfiis _ -fivap Plal/3 /, Q \ 

PL PL pizvr 



The decomposition in Eq. (18) covers all possible angle dependences which occur in the 
process with a single polarized quark. This decomposition gives rise to the decomposition 
of the differential cross section according to 



= 1(1 + cos 2 9)aP + ^ sin 2 9aP + ^co^a) r ' . 
a cos 9 8 4 1 

3 . „ „ rp\ 3 



+ ^ sin 9 cos 9a { p + - sin 9a { p (20) 
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3 Unpolarized and polarized structure functions 



the mean longitudinal polarization of the quark (i.e. the average over the polar 

Ki|j the polar angle dependence was determined, 



In Ref. [ 

angle 9) was calculated, while in Refs. 
including beam polarisation effects. The polar angle dependence of the transverse polar- 
isation was considered in Ref. ||. In this section I present the results and show a few 
instructive steps of the calculations. 






Figure 1: Born term contributions to the single spin case 



3.1 Born term contributions 

The Born term contributions are easy to calculate. The Feynman diagrams for the hadronic 
part (in case of the single polarized case) are shown in Fig. |I|. The unpolarized contributions 
are given by 

Hlj(born) = 2N c q 2 (l + v 2 ), H l L (born) = N c q 2 (l - v 2 ) = H 2 L (born), 
Hl{born) = 2N c q 2 {\ - v 2 ), H^(born) = 4N c q 2 v. (21) 

where N c is the number of colours. The longitudinally polarized contributions read 

Hff(born) = AN c q 2 v, H^(born) = 2N c q 2 (l + v 2 ), 

Hfiborn) = 0, H 2 F e {born) = 2N c q 2 {l - v 2 ). (22) 

For the transverse polarization we finally obtain 

Hj L (born) = -2N c q 2 vJl, H^^born) = -2N c q 2 ^ = H 2 /(born), 

Hl N (born) = 2N c q 2 Jl (23) 

where £ = Am 2 /q 2 = 1 — v 2 and "_L" and "iV" indicate the transverse polarisations per- 
pendicular and normal to the event plane. 
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3.2 First order loop contributions 

The calculation of loop contributions can be formally be done by using the vertex form 
factors, i.e. by replacing the vector and axial vector vertices according to 



r 



1 + A)-f + B{ Pl - p 2 r, 7^75 - ((1 + C)^ + D{ Pl + p 2 y) 75 . 



(24) 



The squared matrix element including first order loop corrections contains the Born term 
contributions as well as first and second order contributions. Only the first order con- 
tributions, as shown in Fig. |2|, are used for a first order calculation. The non-vanishing 
contributions read 



Hi 



H 



it 



2( 



H 



H 1 / 



H} N 
Hf 
Hi N 



oop) = AN c q 2 ({l + v 2 )ReA-2v 2 ReB) , 

oop) = AN c q 2 ((l-v 2 )ReA + 2v 2 ReB) , 

oop) = 2N c q 2 ((1 — v 2 ) ReA + v 2 Re I?) = H 2 L (loop), 

oop) = 8N c q 2 v (Re A — ReB) , 

oop) = 8N c q 2 v (Re A — ReB) , H A L \loop) = 0, 

oop) = AN c q 2 ({l + v 2 )ReA-2v 2 ReB) , 

oop) = AN c q 2 ((l-v 2 )ReA + 2v 2 ReB) , 

oop) = -2N c q 2 v^(l + Z)ImB/Z, 

oop) = -2N c q 2 v^(2ReA+(l-30ReB/£), 

oop) = -2N c q 2 ^l{ReA-v 2 ReB/() = H 2 /(loop), 

oop) = -2N c q 2 vJi(l-i)lmB/i = H 2N (loop), 

oop) = 2N c q 2 v^(2ReA+ (1 - 3£) Re5/£) , 

oop) = 2N c q 2 v^(l + OlmB/^ 



(25) 



the vertex form factors A and B are given by 



Re A 



a s CF 




2 + 




+ 4 



v 



ReB 



a s C F 1 -v 2 



(26) 



47T V 



were we already have replaced the axial- vector form factor C by A according to C = A—2B, 
and where m 2 G = Aq 2 is the squared mass of the gluon which we use for the regularization 
of the first order tree contributions and which is replaces the parameter e = (4 — D)/2 of 
the dimensional regularization by means of 1/e — + ln(47r/i 2 /m 2 ) — > ln(Aq 2 /m 2 ). 



3.3 First order tree graph contributions 

Even though ending up with a different final state, namely a three particle state including 
quark, antiquark and gluon, the first order tree graph contributions are needed for a full 
first order result. For vanishing gluon momentum the (soft) gluon cannot be resolved 
by the detector, the event is counted as two particle final state. In this soft region the 
infrared singularity of the three graph contribution cancels the infrared singularity of the 
loop contribution. This is a general result formulated in the Lee-Nauenberg theorem ||. 

The hadronic part of the diagrams is shown in Fig. ^. The additional degree of freedom 
given by the gluon leads to two additional integration measures, expressible in the energy 
type variables y and z mentioned earlier. The three particle phase space for some specified 
center-of-mass energy y/q 1 is shown in Fig. Integrals to be calculated are of the general 
type (p and q take integer values in a finite range around zero) 




Figure 3: tree graph contributions to the single spin case 




(27) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

y 

Figure 4: Example for the phase space for the three particle decay 



where the integration limits are given by the phase space boundary, 

y- = yfot + K v+ = l-yfi (28) 
z ± (v) = ^{i-y-^ + A + ^i^-A) 2 -^ 1 -^) 2 -^} ( 29 ) 

where VAq 2 is again the gluon mass. While the integration over z can be performed 
easily, the second integration is not so easy to perform, the most complicated (logarithmic) 
integrals are of the shape 




(30) 
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Two issues should be mentioned here. First, the gluon mass is necessary only for regular- 
izing the integrals. Therefore, the parameter A need to appear only at places where the 
integral at the lower limit would lead to singular expressions if the parameter were absent. 
A lot of integrals, therefore, can be calculated for A = 0. On the other hand, the singular 
integrals can be calculated by subtracting and adding a (simpler) integral which has the 
same singular structure. According to 

/'(£, A) = (/'(£, A) - 7^, A)) + A) = A/'(£, A) + A), (31) 

for the difference term A = can be chosen. The integral I' Q can be obtained from I' 
by taking approximations for the integrand close to the lower (infrared singular) limit, 
y = VM + A. So 

*t -> 7 ((2 - Ov ± ^(y - A) 2 - a^v/w) (32) 

and a substitution y — A + y^^T cosh ( helps for the integration. For the integrals with 
A = 0, on the other hand, a combined substitution y — 1 — \/£(l + z 2 )/(l — z 2 ) works to 
separate the integral into integrable pieces. The results can be finally expressed in a very 
compact form by using the so-called decay rate terms. At this point I only show the results 
for the unpolarized case for reasons of brevity, the spin-dependent parts are only different 
but not more complicated, (N = a s N c Cj /4ttv) 

Hh(a s ) = N{2{2 + 7£)v + (48 - 48£ + 7£ 2 )* 3 + 2^ (2 - 7f + ^ (2 + 3f)) * 4 

-2£(2 + 30*5 - 4(2 - ((2 - £)(*8 - h) + 2v(t 10 + 2t 12 )) }, 
Hlias) = ^{12^ + 2(6-0*3 + 2^(1-^4 

+2^*5 - 4 ((2 - 0(ts - * 9 ) + 2u(*io + 2t 12 )) }, 
Hl(a s ) = N{~(W - 46£ + 3e)v + ^(88 - 32£ + 3£ 2 )t 3 

-2^ (2-7^ + y^(2 + 30) t 4 + 202 + 30*5 

-2£ ((2 - 0(^8 - t 9 ) + 2v(t 10 + 2t 12 )) }, 
H 2 L (a s ) = ^{^(10-0^ + j(24-16e-3f)t 3 

-2^(1 - v£)*4 - 2^5 - 2 ((2 - 0(*s - t 9 ) + 2v(t 10 + 2t 12 )) }, 
ffj.(a.) = -8fiVwr, 

Hf,(a.) = iV{-16 v /e(l- V / e)-16(ti-t 2 ) + 8(2 -30^3 

-4(4 - 50*6 - 8v ((2 - 0(*s - t 7 ) + 2u(tio + tn)) } (33) 

where the decay rate terms, containing logarithms and dilogarithms, are found in Ref. fllOl 
together with the polarized contributions. More interesting at this point are the results we 
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obtain. In Fig. ^ (top) the total cross section for pair production of top quarks is shown in 
dependence on the center-of-mass energy y/q 2 ', with and without the first order radiative 
correction. The polar angle dependence of the longitudinal polarization 

= <Mt) -dcr(l) 

do{]) + da{l) { ' 

of the top quark is shown in Fig. |6] (bottom). In Fig. [7| the polar angle dependence of 
the normal (top) and perpendicular polarization (bottom) of the top quark is displayed. 
What is not shown here is the change with respect to the Born term contribution. It 
figures out that despite the fact of rather large corrections to the total cross section, the 
corrections to the polarizations are quite small (about 5%). A detailed discussion of the 
results can be found in Refs. || P, pLQ| . Quantities like the polarization can be measured 
by the asymmetry of the decay products of the top quark. 

3.4 Finite gluon energy cut 

Besides the comparison with the massless case, resulting in anomalous spin-flip contri- 
butions which will be shown later, additional work has been investigated in the exact 
calculation of the contributions for a given gluon energy cut, dividing the phasespace up 



into a soft and a hard part [16]. This "excurse" will not be taken here. 



4 Longitudinal spin correlation 

If both quark and antiquark are considered to be polarized, the correlation between their 
spins can be measured. In Ref. [jTB) the correlation of longitudinal polarized quark and 



antiquark is considered, in Ref. |L4] the angular dependence is considered as well. Because 
of the representation 



^(V(i-y) 2 - £;o,o,i- y ) (35) 



i 



'AT = -^=(\/(l-z)i-Z;(l-z)sm6 12 ,0,(l-z)cos6 12 ) 



for the longitudinal spins where 



yz + y + z-l + £ 
cos #i2 = , , , (36) 

:i-y) 2 -W(i-*) 2 -£ 



there occurs an additional square root y (1 — z) 2 — £ in the denominator which makes the 
integrations more complicated. However, the same integration methods work also in this 
case, an appropiate substitution could be found to make the integration feasible. New 
decay rate terms had to be defined in order to describe the final result (see Ref. |b 
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4.1 Result for the correlation 

The correlation is denned as 

r = da(W-da(n)-da(n) + da(U) 

do-(TT) + MU) + da(i^) + da(ii) 1 ' 

The polar angle dependence of this result is shown in Fig. ||. 

4.2 Anomalous spin-flip contributions 

It is an easy task to calculate results for polarized structure functions also in the case 
where the quark mass is assumed to be massless, even though this calculation have to 
be performed in dimensional regularization. The results, however, differ from results one 
obtains by taking the results for the massive quark as obtained earlier and calculating the 
limit m — > 0. By adding in the Born term contributions one finally obtains in the m — > 
limit {Cp = 4/3 is made explicit here) 



4 a s t t 

o X S 1 S 2 
3 7T 



4 

\ 3 7T 
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H F [s 



1) °2) 



H, 



tA d l> b 2) 



'1> 



N c q\s[ - 
N c q 2 
N c q 2 





[2 a." 




fl + Ox^- 


- x — 


) 


V 7T 


.3 71 _ 





Si - s. 



1 , 1 a * 
1 + - X 

3 7T 



2 a s 
- x — 

3 7T 



Si - Si 



+ |x^ + [0] 

O 7T 



7T 



4 

- X 
3 7T 



s £ £ 

s 1 s 2 



(3f 



where the anomalous spin-flip contributions are made explicit by the square bracket no- 
tation. These contributions have their origin in the collinear limit where the spin-flip 
contribution proportional to m survives since it is multiplied by the 1/m collinear mass 
singularity. Because the anomalous spin-flip terms are associated with the collinear sin- 
gularity, the flip contributions are universal and factorize into the Born term contribution 
and an universal spin-flip bremsstrahlung function [|17]]. This explains why there is no 
anomalous contribution to and why the anomalous flip contributions to Hy and H F , 
for instance, are equal. In fact the strength of the anomalous spin-flip contribution can 
directly be calculated from the universal helicity-flip bremsstrahlung function listed in 
Ref. [O. 



5 Gluon polarization 



The polarization of gluons in e + e~-annihilation [18|, |l9j, in deep inelastic scattering p0 | 
and in quarkonium decays |TR ETJ has been studied in a series of papers dating back to 
the early 80 's. Several proposals have been made to measure the polarization of the gluon 
among which is the analytic proposal to angular correlation effects in the splitting process 
of a polarized gluon into a pair of gluons or quarks. 

Results for the linear and circular polarization of gluons produced in the e + e~ anni- 
hilation process for massive quarks with subsequent gluon emission are presented in this 
section. Similar to the joined quark- ant iquark density matrix the two-by-two differential 
density matrix d<sr = d<j\ G \' G of the gluon is used with gluon helicities Ac = ±1 in terms of 
its components along the unit matrix and the three Pauli matrices. Accordingly one has 

de = -(dat + da x v x + do y v y + da z sr z ), (39) 
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where da is the unpolarized differential rate and da = (da x , da y , da z ) are the three com- 
ponents of the (unnormalized) differential Stokes vector. After azimuthal averaging the 
^/-component of the Stokes vektor da y drops out, one is left with the x- and ^-components 
of the Stokes vector which are referred to as the gluon's linear polarization in the event 
plane and the circular polarization of the gluon, respectively. 

5.1 The differential cross section 

The differential unpolarized and polarized cross section, differential with regard to the 
polar beam-event orientation and the two energy-type variables x = 2p^ ■ q/q 2 and w = 
2(pi — P2) ■ q/q 2 (with q = p 1 + p 2 + P3) are then given by 

drr^ 3 / da 1{x) dcr 2 ^ 

:^^9)[ gil p-+g l2 p-\ , (40) 



d cos 9 dx dw 8 \ dx dw dx dw 

3 . , J dal (x) da 2 L {x) \ 3 „ daf ] 
+ -sm 2 # \g u — L — + gi2-r L -r- + 7 cos9 g u -^—, 

4 I dx dw dx dw J 4 dx dw 

da z 3^ 2n , daff 3 „ / dal z daf \ 

l + cos 2 e)g 1A -^ + -cose\g Al —^- + g A2 —^-\. (41) 



d cos 9 dx dw 8 dx dw 4 y dx dw dx dw ) 

The notation da^ stand for either da or da x , and the same for da l ^ x \ Of course, in this 
case no Born terms occur, the first order tree diagram is the leading order contribution. 
The phase space in terms of the variables x and w is similar to the one in y and z (see 
Fig. U), standing on the sharp corner. The phase space limits are therefore symmetric, 

X- = 0, x + = 1 — £, 



w ± (x) = ^J - (42) 

5.2 Results for the gluon polarization 

The integration can be done also in this case by using an appropiate substitution (see 



Refs. [[11], 0). However, it is also worth to consider the results still depending on the 



gluon energy variable x. The results are given by 

4 1 2 

H^jix) = -32(--4 + x)-w + (x)+32(--2 + a;)4(x), 

\ nf / nr \ nr / 

Hl(x) = l6(4 + e)(--l)-w+(^)-8e(— -2-a;)4(a;), (43) 
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where 



H 2 L %x 



H 4 F (x) = 64(— ^-l) -32(^-2 + x)t e _(x 
Hj?(x) = -128(--l)-w + (x) + 64:(--l)t e + 



Hff(x) = -32e(--l)-«; + (^) + 16e(— -2)ti(x), 

00 00 00 

Hf{x) = 16(4 + (- - l) - 8^(^-2)4^), (44) 

= -iee(i-i)W)+ se(^-2)4(x), 



Hp(x) = 64(i-i-l)-64(l-l)tl(x), 



;r 

H%{x) = 64(1 -£ - x) -32(2- x)tl(x), #| z (x) = 0, 

M a (x) = -32(4 -3x)-w+(x) +32(2- x)ti(x), = (45) 

x 



»+<*)= *V^r- ^» =lm U-»-vi-»-J ' ti(l)=h, (— J- 

(46) 

For different values of x in terms of x max = 1 — £ the polar angle dependence of the linear 
and circular polarization of the top quark is shown in Fig. R 



6 Conclusion 

First order radiative QCD corrections for polarization observables can be done analytically 
and are done in most of the cases. The corrections for the polarization observables figure 
out to be small. This, however, proves the quality of perturbation theory estimates. The 
polarization observables can be measured by the asymmetry in the decay rate of final state 
leptons. The consideration of massive quarks are important especially at the high energies 
which are reached by the NLC. The calculation techniques can be combined with other 
techniques like the high spin tools developed in the theory group at Tartu University, a 
future collaboration is scheduled. 
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Figure 6: total cross section dependence on the center-of-mass energy y/q 1 with and without 
first order radiative corrections for the top quark pair production (top) and polar angle 
dependence of the longitudinal polarization for the center-of-mass energies 380 GeV (solid), 
500 GeV (dashed) and 1000 GeV (dashed-dotted) including radiative corrections 
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Figure 7: angular dependence of the normal polarization (top) and the perpendicular 
polarization of the top quark (bottom) for the center-of-mass energies 380 GeV (solid), 
500 GeV (dashed) and 1000 GeV (dashed-dotted) including first order radiative corrections 
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Figure 8: polar angle dependence of the linear (top) and circular polarization (bottom) of 
the gluon in the case of top quark pair production 
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